x <- c(5)
y <- c(3)
alpha <- 4E-2
beta <- 4E-3
gamma <- 4E-3
delta <-4E-3


for (i in 1:4000) {
  lastx <- x[length(x)]
  lasty <- y[length(y)]
  dx = lastx*(alpha - beta*lasty)
  dy = -lasty*(gamma - delta*lastx)
  x <- c(x, lastx + dx)
  y <- c(y, lasty + dy)
}

plot(x, type="l", col="green",ylim=c(0,max(x)), xlab="Dienas", ylab="Zveri")
lines(y, type="l", col="black")

#plot(x,y, type="l", col="gray", xlab="pavi\u0101ni",ylab="gepardi")